The European Physical Journal Plus

Micro-elastography is an optical technique that studies elastic waves for the mechanical characterisation of micrometric objects, such as cells. We propose to adapt this technique for the characterisation of millimetre-sized samples using a white light microscope. The objective is to perform a rapid, global characterisation of the elasticity of a biopsy. The millimetre-sized samples to be characterized are embedded in an agarose gel. A vibrator generates shear waves in this gel that transmit naturally inside the sample. This technique removes the need for precise manipulation of the wave source. A high-speed camera records the propagation of the waves in the sample. Their velocity is calculated using a noise correlation approach. Due to the lack of millimetric phantoms of calibrated elasticity, we choose to validate this method with a three step process. The experimental setup is first validated on homogeneous gels, then on biological samples of increasing elasticity, biopsies of beef liver hardened by heating, and finally on biological samples of clinical interest: biopsies of mouse endometrium. This method can be applied to all types of biological tissue, paving the way for rapid mechanical characterization of biopsies.

During embryogenesis, as cells proliferate and assemble into tissues, they undergo a sequence of transitions between distinct molecular states eventually giving rise to a cellular population consisting of an appropriate distribution of specific functional cell types. Recent progress on the dynamics underlying decision-making in developmental landscape makes it feasible to start analysing the amount of information involved in constructing such systems. To explore this we introduce the notion of potency of a developmental landscape and attempt to calculate it for two development systems of current interest, in-vitro differentiation of epiblast-like cells into neural and mesodermal progenitors and the worm vulva patterning system. Our approach integrates concepts from developmental biology, information theory and dynamical systems to estimate both the number and identity of signalling regimes that give rise to distinguishable temporal response patterns.

The light absorption of the protein cryptochrome and its chromophore FAD is important for the regulation of circadian rhythms and in some species for sensing magnetic fields. To compute the absorption spectrum of chromophore, typically only a small region is treated quantum-mechanically due the high computational cost of spectroscopic calculations. We present a formalism that allows a quantum-mechanical treatment of not only the chromophore but also the neighbouring amino acids which differ from species to species. This is achieved by using the real-time time-dependent Hartree-Fock method. This method allows extending the quantum domain from typically only a few dozen atoms up to around 1,200 atoms for the largest calculations. The presented framework allows the treatment of neighbouring tryptophan residues or the cofactor molecule MTHF in the same calculation and allows to extract information of which regions absorb light depending on wavelength. The presented results also show that the environment around the chromophore FAD amplifies the light absorption in cryptochrome.

Vertebrate gastrulating mesoderm is a prototypic example of a mesenchymal-like tissue undergoing extensive remodelling. While the tissue may be globally represented as a viscoelastic material, the actual biological material is intrinsically complex. To get to a real understanding of its properties, one needs to move to the mesoscale, linking cellular properties to collective phenomena. Vertebrate embryos also display a remarkable variability in mechanical properties, despite which they robustly complete gastrulation. This study attempts to explore these aspects by dissecting Xenopus mesoderm cell behaviour in a minimal system, using aspiration through a microfluidic system to impose controlled stress to a mesoderm aggregate. We show that beyond estimating global rheology at the tissue scale, it is possible to infer a wealth of information based on cell morphology and dynamics. Our data are consistent with collective behaviour being mostly dictated by the balance between the capacity of cells to stretch and the resistance to cell-cell contacts, which limits cell-cell intercalation and thus tissue remodelling. Importantly, tissues are not only able to transmit stress over a distance, they also clearly react to it through actively reinforcing cell-cell mechanical coupling. This adaptative property is found through a broad range of tissue stiffness, and adhesion strength appears to scale with the elastic modulus, suggesting that cell stiffness may ultimately be the key parameter setting mesoderm rheology and accounting for the large differences observed between embryo batches.

Ultraweak photon emission is the spontaneous emission of extremely low levels of light from a broad range of biological systems. Recent studies have reported that UPE measured extracranially can serve as a potential non-invasive biomarker of brain activity. Here, we show that this interpretation suffers from serious problems. First, when observed under properly dark conditions, the UPE from the head is much weaker than what is reported in certain papers on brain UPE from human heads. Signals detected in these studies are overwhelmingly dominated by background light. Second, photons at wavelengths < 600 nm are strongly attenuated by scalp and skull tissues, and longer wavelengths fall largely outside the effective spectral sensitivity of the photomultiplier tubes (PMTs) used. As a consequence, even if UPE from the head is detected under properly background-free conditions, it is likely to be dominated by emission from the scalp rather than from the brain, certainly as long as PMTs are used. Our results emphasize the importance of careful experimental design to make genuine progress on this important question.

The Poisson-Nernst-Planck (PNP) system is an accurate model of electrodiffusion of ionic species. It is commonly used in situations where nanoscale resolution is required, for instance close to ion channels in the membranes of biological cells. The inherent stiffness of the equations has made them challenging to solve and has limited the applicability of the system. In particular, the time step required for stable solutions has typically needed to be very short (nanoseconds), which makes simulations on the time scale of an action potential (milliseconds) difficult. Recently, it has been observed that avoiding operator splitting and instead solving the concentration equations and the electrostatic equation in a coupled manner relaxes the time-step limitation considerably. However, no theoretical explanation of this observation has been provided. Here, we aim to explain why the coupled scheme allows much larger time steps. We illustrate the mechanism by considering special cases that define necessary, but not sufficient, conditions for stability. We also show that these conditions remain relevant for the fully coupled PNP model in 3D.

This work focuses on the global and partial identification problem for fractional differential equations. We provide a general numerical procedure based on global and local optimization algorithms with two refinements for biological systems that ensure solution positivity and homogeneous parameter units. The method is applied to a new fractional model of Dengue outbreak called the Fractional Homogeneous Nishiura (FHN) model, calibrated using data of newly infected people in Cape Verde. We show that our identification method yields a better fit between data and model solutions than previous approaches and that our FHN model captures the dynamics of Dengue more closely than existing systems.

A major challenge in neuroscience is to predict how currents in nanodomains affect voltage and ionic concentrations. Cable and Rall theory provide analytic current-voltage relations by neglecting concentration gradients, and the impact of concentration gradients is usually studied numerically with the Poisson-Nernst-Planck (PNP) model. A precise quantitative understanding of the combined dynamics remains limited because analytic current-voltage-concentration relations are missing. In this work we derive such relations using a novel approach based on cross-diffusion equations. For narrow cylindrical domains, we derive time-dependent and steady-state expressions that explicitly show how currents affect voltage and ionic concentrations. We find that the influx of only one ion can significantly change the concentrations of all the other ions even if no channels for these ions are present. After a current injection we compute a biphasic voltage transient where the small-time asymptotic corresponds to the steady-state solution of the cable equation. We show that the accuracy of cable theory prediction for the voltage depends on how the current is distributed among the various ions. Finally, we develop an iterative method to accurately compute steady-state profiles for voltage and concentrations using first-order results by subdividing a cylinder into small segments.

Biological tissues grow at rates that depend on the geometry of the supporting tissue substrate. In this study, we present a novel discrete mathematical model for simulating biological tissue growth in a range of geometries. The discrete model is deterministic and tracks the evolution of the tissue interface by representing it as a chain of individual cells that interact mechanically and simultaneously generate new tissue material. To describe the collective behaviour of cells, we derive a continuum limit description of the discrete model leading to a reaction-diffusion partial differential equation governing the evolution of cell density along the evolving interface. In the continuum limit, the mechanical properties of discrete cells are directly linked to their collective diffusivity, and spatial constraints introduce curvature dependence that is not explicitly incorporated in the discrete model. Numerical simulations of both the discrete and continuum models reproduce the smoothing behaviour observed experimentally with minimal discrepancies between the models. The discrete model offers further individual-level details, including cell trajectory data, for any restoring force law and initial geometry. Where applicable, we discuss how the discrete model and its continuum description can be used to interpret existing experimental observations.

Ultraviolet (UV) radiation induces DNA damage associated with cancer and aging, yet the sequence dependence of UV absorption remains to be investigated. Here, we present a systematic study of the UV absorption spectra of DNA based on all-electron Hartree-Fock calculations. We analyze all possible sequences up to four base pairs, as well as longer randomized sequences and genomic nullomers - motifs which are missing in a given genome. We observe a pronounced sequence dependence: cytosine- and guanine-rich motifs exhibit significantly enhanced absorption, whereas adenine-thymine-rich sequences absorb up to four times less in the mid-UV range. Notably, the human genome is biased toward adenine-thymine-rich sequences, giving it an increased susceptibility to UV-induced damage. In addition, we introduce a computational framework enabling spectral calculations of large DNA and RNA fragments, opening the door to large-scale optical analyses.

Airway epithelium plays a major role as the primary interface between human body and the external environment, acting both as a physical and functional barrier. In vitro airway models that reproduce the epithelium architecture are therefore a valuable tool for studying infection, inflammation, and transport processes. In this work, we present a label-free, non-invasive method to visualize and measure mucociliary transport in air-liquid human models using third-harmonic generation (THG) microscopy with an optical parametric amplifier laser source at 1300 nm. By exploiting the intrinsic nonlinear contrast at optical heterogeneities, THG provides high-resolution images of both epithelial structures and of the overlying mucus layer without the need for fluorescence staining or sample processing. Time-lapse THG imaging reveals depth-dependent transport dynamics within the mucus, offering new insights into mucociliary transport mechanism. Our approach offers a physiologically relevant way to assess mucociliary function in vitro and could support studies on respiratory diseases, drug delivery and efficacy, and epithelial remodeling. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=117 SRC="FIGDIR/small/717621v1_ufig1.gif" ALT="Figure 1"> View larger version (52K): org.highwire.dtl.DTLVardef@62e8acorg.highwire.dtl.DTLVardef@199a8b7org.highwire.dtl.DTLVardef@113bb84org.highwire.dtl.DTLVardef@7be3f8_HPS_FORMAT_FIGEXP M_FIG For Table of Contents Only C_FIG

Transparent and composite surfaces pose a fundamental challenge for stereo photogrammetry: optically smooth glass produces no detectable surface features under visible illumination, making three-dimensional reconstruction impossible without surface preparation. This excludes optical components such as lenses and cover glasses, composite assemblies, and semi-translucent biological specimens from non-contact geometric measurement. Here we show that coherent speckle illumination at 266 nm overcomes this limitation by exploiting wavelength-dependent scatter enhancement, generating sufficient backscattered signal on surfaces that are entirely invisible under visible illumination. We developed a multispectral stereo system and evaluated three illumination modalities under identical acquisition conditions. On transparent glass, both visible modalities produce complete reconstruction failure, recovering only non-transparent holder structures. Ultraviolet speckle illumination at 266 nm enables dense reconstruction of the same surfaces. We demonstrate recovery of an uncoated plano-convex lens with a fitted radius of 30.946 mm and point-cloud standard deviation of 106.5 {micro}m, defect detection on a transparent cover glass without surface preparation, and reconstruction of a semi-translucent biological specimen. On metrology-grade reference objects, ultraviolet speckle achieves a standard deviation of 116 {micro}m and completeness exceeding 93%, approaching the performance of optimised visible structured illumination. These results establish ultraviolet speckle photogrammetry as an enabling approach of optical metrology to otherwise uncooperative surfaces, with relevance to optical manufacturing inspection and biological surface analysis.

With the development of microfluidics, it has now become possible to assess the susceptibility of bacteria to antibiotics at the single-cell level instead of relying on population measurements. Such studies are particularly relevant when the growth of bacterial population in the presence of antibiotics is heterogeneous. Here, we build a model to describe such a case, and apply it to experimental measurements on a small population of E. Coli exposed to ciprofloxacin, a drug which is well known for triggering a bistable response.

In tissue engineering, it is important to conceive and construct artificial bio-mimetic scaffolds able to foster cell migration as this is a fundamental process in wound healing and tissue regeneration. In order to do that, cubically symmetric and triply periodic porous structures have been identified as promising candidates for instance for the reconstruction of artificial cartilages and bones, also due to their tunable mechanical characteristics and highly inter-connected porous architectures that mimic the trabecular bone hyperboloidal topography. We propose here a mathematical approach that might be helpful to identify what are the best geometrical characteristics of such scaffolds, in order to promote cell migration into the porous structures and speed-up their re-population. The method is based on the observation that cell nucleus deformations should be avoided, yet assuring a good possibility for the cell to reach the wall of the porous structure. Mathematically speaking, this leads to the problem of identifying the size of the largest sphere that can pass, without being stuck, through the pores of the bio-mimetic scaffold.

Mass-conserved reaction-diffusion systems are used as mathematical models for various phenomena such as cell polarity. Numerical simulations of this system present transient dynamics in which multiple stripe patterns converge to spatially monotonic patterns. Previous studies indicated that the transient dynamics are driven by a mass conservation law and by variations in the amount of substance contained in each pattern, which we refer to as "pattern flux". However, it is challenging to mathematically investigate these pattern dynamics. In this study, we introduce a reaction-diffusion compartment model to investigate the pattern dynamics in view of the conservation law and the pattern flux. This model is defined on multiple intervals (compartments), and diffusive couplings are imposed on each boundary of the compartments. Corresponding to the transient dynamics in the original system, we consider the dynamics around stripe patterns in the compartment model. We derive ordinary differential equations describing the pattern dynamics of the compartment model and analyze the existence and stability of equilibria for the reduced ODE with respect to the boundary parameters. For a specific parameter setting, we obtained results consistent with previous studies. Moreover, we present that the stripe patterns in the compartment model are potentially stabilized by changing the parameter, which is not observed in the original system. We expect that the methodology developed in this paper is extendable to various directions, such as membrane-induced pattern control.

Finite element (FE) head models are valuable tools for investigating brain injury mechanics, with their reliability critically dependent on accurate material modelling. White matter (WM) is often considered mechanically anisotropic due to its aligned axonal fiber architecture and is commonly represented using fiber-reinforced hyperelastic formulations such as the Gasser-Ogden-Holzapfel (GOH) model. A fundamental assumption of the GOH model is that fibers contribute only in tension and not in compression, requiring the use of tension-compression switches. However, inconsistencies were noted in the formulation of tension-compression switches with the influence on computational biomechanics unknown. To address this knowledge gap, three commonly used switching schemes - differing in both the switching parameter and the treatment of compressed fibers - were theoretically elaborated and numerical implementation within the GOH framework to simulate the mechanical anisotropy of WM in impact simulations. Results from the case-based and group-level analyses demonstrated that both the switching parameter and the treatment of compressed fibers affected WM deformation. Significant cross-scheme strain differences were noted in the first principal strain at the element level and fiber strain at the fiber level. These findings highlighted the mechanical role of tension-compression switch in the GOH-based brain modelling and advocated the adoption of fiber stretch itself as the switching parameter to discriminate the tensile and compressive fibers. The current study provides important guidance for the anisotropic constitutive models in brain tissue and calls for direct verification of the tension-compression switch hypothesis in axonal fibers.

Digital holographic microscopy systems in a common-path configuration, compared to systems with a separate reference arm, offer a compact design and resistance to disturbances. They can operate with partially coherent illumination, reducing speckle noise. However, they are limited by the overlapping of the object beam and its laterally shifted replica. As a result, images from different regions of the object overlap on the detector, preventing imaging of dense samples. We present the wavelength-scanning replica-removal method, which solves this problem by enabling the separation of information from both replicas and thereby doubling the effective field of view (FOV). The wavelength-scanning multi-shear replica removal algorithm plays a key role in reconstructing the undisturbed phase from a series of holograms recorded with variable shears. The shear value is controlled by changing the illumination wavelength. This enabled the development of two measurement modes: time-domain wavelength scanning for high-quality imaging, and a single-shot mode with frame division into color channels to improve temporal resolution. The method was validated using resolution tests and biological samples - neurons and dynamic yeast cultures. By combining the advantages of the common-path configuration with dense-structure imaging and dynamic processes, the proposed method constitutes a versatile tool for quantitative phase microscopy.

Invertebrate vision relies on bistable visual pigments flipping upon photon absorption between rhodopsin and metarhodopsin states. In living butterflies, the UV-VIS absorption spectra of rhodopsin and metarhodopsin, respectively with 11-cis and all-trans isomers of 3-hydroxy-retinal (A3) chromophore, can be conveniently recorded from the eyeshine, the light reflected from the compound eye after passing twice through the light-guiding rhabdoms. * Here, a microscope coupled with a broadband LED source and a microspectrometer was used to record photorelaxations reported in eyeshine reflection spectra. Fitting temporal exponential relaxations to log-reflectance arrays yielded transient and baseline spectra that are analogous to absorbance difference and sum, respectively. Both types of spectra were subjected to singular value decomposition and to fitting of templated visual pigment absorption spectra. * The compound eye of the high brown fritillary Fabriciana adippe was exposed to a series of second-long broadband light pulses, causing photorelaxations with time constants between 40 and 120 ms that led to 80% metarhodopsin in equilibrium. The transient and baseline spectra were fitted with pigment templates, estimating the alpha peak wavelength 547-552 nm for rhodopsin and 496-501 nm for metarhodopsin. The metarhodopsin to rhodopsin alpha peak absorbance ratio 1.25-1.35 is consistent with the isosbestic wavelength at 530 nm. The second isosbestic wavelength indicates that rhodopsin beta (UV) peak absorbs more strongly than metarhodopsin below 405 nm. * Baseline spectra, which were not explicitly analysed in previous studies, enable concatenation of exposures, monitor long-term changes of pigment, and enhance the estimation of beta peak parameters. * The method can be directly used in many butterflies and could be adapted to other insects, particularly fruitflies, facilitating studies of the relation between the visual pigment spectra and the opsin sequences. Spectroscopic results can be complemented with physiologically measured photoreceptor spectral sensitivity datasets and analysed with the same global fitting procedure.

Cerebrospinal fluid (CSF) is a Newtonian fluid that bathes the brain and spinal cord and oscillates in response to the physiological periodic changes in brain volume, of which the cardiac cycle is a major driver. Understanding this motion is essential for clarifying its contribution to solute transport, waste clearance, and drug delivery. In this work, we study oscillatory and steady streaming flow in the cranial subarachnoid space using a lubrication-based theoretical framework. The model represents the cranial CSF compartment as a thin fluid layer bounded internally by the brain surface and externally by the dura, driven by time-dependent brain surface displacements. We first derive simplified governing equations for flow over an arbitrary smooth sphere-like brain surface and obtain analytical solutions for an idealised spherical geometry with uniform displacements. We then incorporate realistic displacement fields reconstructed from MRI measurements in healthy subjects and solve the reduced equations numerically. The results show that oscillatory forcing produces a steady streaming component that may enhance solute transport compared with diffusion alone. This work provides a mechanistic description of the flow generated by physiological brain motion and highlights the potential presence of steady streaming in cranial subarachnoid fluid dynamics.

The spatial resolution of optical imaging systems is fundamentally restricted by the diffraction limit. However, in widefield live-cell microscopy, the achievable resolution is further constrained by the specimen motion, which indicates the existence of a fundamental spatio-temporal resolution trade-off between signal accumulation during the full frame integration and the resulting motion blur. To improve the fidelity with which moving objects can be imaged, a quantitative understanding of this spatio-temporal trade-off is necessary. Here, we present a systematic analysis of motion-induced resolution dynamics measured with spectral signal-to-noise ratio (SSNR). We developed a simulation framework which models the image formation of objects undergoing arbitrary motion, to evaluate the degradation of the spatial resolution under translational and rotational dynamics. Our results demonstrate that for translating objects, the spatial resolution is anisotropically reduced as a function of the orientation of the object relative to the motion vector, leading to the spectral signal-to-noise ratio degrading by up to 50% and the resolution by up to 40% for a 90{degrees} change in the motion direction. Furthermore, we show that for rotational motion, conventional radially averaged metrics such as the Fourier Ring Correlation are not able to quantify the effects of angular blur. On the other hand, the SSNR is able to accurately quantify this degradation. These findings underscore the necessity of an object-oriented imaging approach, in which acquisition parameters such as exposure time are tuned to specific biological spatio-temporal characteristics to optimize the trade-off between motion blur and spatial fidelity.